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Abstract 

A study of the reaction tt + + d — > p + p has been performed in the energy 

range of 18 - 44 MeV. Total cross sections and differential cross sections at 

six angles have been measured at 15 energies with an energy increment of 

1-2 MeV. This is the most systematic data set in this energy range. No 

structure in the energy dependence of the cross section has been observed 

within the accuracy of this experiment. 
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I. INTRODUCTION 



Total and differential cross sections for the reaction ir + d — > pp have been measured 
with high accuracy at pion energies above 100 MeV. Currently, experimental efforts in this 
energy range are aimed at measuring spin dependent observables. The situation at pion 
energies below 100 MeV is less complete, especially for energies of a few tens of MeV [|TJ. 
The disagreement between the experimental data for the total cross section in some cases 
exceeds 20%, much larger than the quoted uncertainties. The experimental data show a 
deviation from smooth behavior in the total cross section of the reaction 7c + d — > pp close 
to T n ~ 30 MeV (y/s = 2.04 GeV/c 2 ) f2|. More recent measurements of the partial total 
cross section for the inverse reaction at SATURNE have verified this effect ||. In a recent 
experiment of the Dubna-Gatchina group on pion absorption on carbon Q a dip in the 
energy dependence of the quasi-deuteron component of absorption near T n =28 MeV has 
been observed. The measurements of pion absorption on the deuteron || performed by the 
same group also indicated possible structure at a pion energy of 30 MeV. A more recent 
experiment at LAMPF || did not observe any dip on carbon, but these data suggest a 
change in the angular distribution of protons at pion energies near 25 MeV. Some possible 
alternative explanations of the structure in this energy range are given in Ref. [0,f|]. One 
of them attributed the structure to the excitation of a diproton resonance in the 3 P2 NN 
state. However, these experiments are not accurate enough to reach a conclusion concerning 
the structure under the discussion. A satisfactory solution fitting the available database has 
not been found in partial wave analysis, suggesting underlying systematic problems with the 
data base. 

The experiment described here was intended to resolve the discrepancies in existing data 
and to verify the observed structure in the excitation function of pion absorption on the 
deuteron at low pion energies. The differential and total cross sections of the n + d — > pp 
reaction have been measured with fine steps in incident pion energy. 
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II. EXPERIMENT 



A. Pion Beam and Target 

The experiment was performed at the Low Energy Pion channel (LEP) of the Clinton 
P. Anderson Meson Physics Facility (LAMPF). Positive pions with energies of 21, 23, 25, 26, 
27, 28, 29, 30, 31, 33, 35, 37, 39, 41 and 45 MeV were used. The pion beam had an average 
intensity of a few 10 4 /sec, allowing the beam particles to be counted while not overloading 
the BGO detectors. The pion fraction varied from 78% for 45 MeV beam to 33% for 21 MeV 
beam. The momentum bite of the LEP channel was set to 1% for the highest energies and 
increased to 4% for the lowest ones in order to increase the pion flux. 

The target was composed of CD 2 with a cross sectional area of 1 x 1cm 2 and an areal 
density of 0.469 g/cm 2 . It was attached to a thin paper pipe and placed in the center of 
the BGO ball. The supporting pipe was aligned along the beam axis. The diameter of the 
pipe was big enough to keep its walls out of the beam. A 0.25-mm-thick plastic scintillator 
SI with a cross section of 6 x 6 mm 2 was located just before the target. Downstream of 
the target a rectangular array of nine Csl scintillators DA1-DA9 was preceded by a 10-mm- 
thick plastic scintillator, S2. DA1-DA9 and S2 could be used for the detection of the most 
forward going reaction products as well as for a determination of beam composition. In this 
experiment only the central detector, DA5, of the array was used. A coincidence between SI 
and the central detector, DA5, of the downstream array was used as a beam monitor. Every 
1000'th beam event (event triggered by a SI ■ DA5 coincidence) was read out for further 
analysis to determine the pion fraction which was needed for absolute normalization of the 
cross sections. 

B. BGO Ball Spectrometer 

A large solid angle detector, the LAMPF BGO ball, was used to detect the reaction 
products in this study. Detailed information on the BGO ball can be found in Ref. [|7||§. The 
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BGO ball consists of 30 phoswich detectors. The detectors of the array were of pentagonal 
and hexagonal shape and tightly packed to form a truncated icosahedron of 32 sides. Two of 
the 32 sides are opened for the beam entry and exit. The detectors were distributed about an 
inner radius of 6.1 cm from the center of the array to the center of each crystal face, and were 
arranged in six groups centered at laboratory scattering angles of 9 = 37°, 63°, 79°, 102°, 116°, 
and 142°. Each detector had a solid angle of about ^ x 4-7T sr and was supported in a 0.5- 
mm-thick electro-formed nickel can which had a 0.05-mm-thick entrance window. Each 
detector consisted a 3-mm-thick NE102 plastic scintillator optically coupled to the front of 
a 5.6-cm-thick bismuth germanate (BGO) crystal, with a 7.62-cm-diameter photomultiplier 
tube on the back. Since the decay constant of the BGO scintillator is much longer than that 
of the plastic scintillator (250 ns vs 1.5 ns), the anode signal was time sliced to provide both 
AE (fast) and E (slow) signals for charged particle identification (pions, protons, deuterons, 
etc.), and for identification of neutrons and gamma rays. The crystals were thick enough to 
stop up to 185-MeV protons and 90-MeV pions. The time resolution of the detectors was 
about 1 ns, sufficient to eliminate events with hits from different beam bursts (the LAMPF 
beam has a 5-ns microstructure) . The light output of BGO scintillator depends significantly 
on the temperature of BGO material |Q. To minimize fluctuations in temperature of the 
BGO, a tent-like structure was built to isolate the BGO ball from its surroundings. 

The event trigger consisted of a coincidence between the target detector, SI, and at least 
one BGO crystal in anti-coincidence with DA5. 

III. DATA ANALYSIS PROCEDURE 

A. Raw Data Processing 

The raw data for each event contains information about the energy deposited in the 
plastic and the BGO for all detectors and timing information with respect to the beam 
counter. In the first step of the analysis a time gate was applied to all of the signals to remove 
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accidentals from further analysis. The next step was the determination and application of 
the constants used for unmixing the AE — E information from the phoswich detectors. 
The signal from each phoswich detector was integrated in two different ADC channels with 
different time gates, 50 and 250 ns long. The AE and E information was separated using: 

AEi = (dEi — Eli ■ Ri,i — ■ DEGAINi , (1) 

Ei = (Eli — dEi ■ R ii2 — Z ij2 ) ■ EGAINi , (2) 

where dEi and Eli represent the raw data from the ADC's with the short and long gates, 
respectively, and DEGAINi, EGAINi are the coefficients for conversion of raw ADC data 
to energy in MeV. The relative fraction of a long signal in a short gate and vice versa are 
given by the mixing parameters, R it i and i? ij2 , respectively. The R+j and offsets, Z it j, were 
determined using a least square fit for the events along AE and E axes, where only one 
component of a signal (short or long) exists. These are produced by charge particles which 
stop in the plastic and leave a AE with no E signal, or by neutral particles which interact 
in the BGO and leave a E with no AE signal. An example of a AE — E distribution from a 
phoswich is shown in Fig. [I]. One can clearly see regions corresponding to protons and pions. 
The outlined region along the E axis is due to neutral particles and the outlined region along 
the AE axis is due to short range charged particles stopped in the plastic scintillator. 

B. Calibration and Stabilization of the BGO Ball 

An initial energy calibration of the BGO ball (DEGAINi and EGAINi) was made by 
using two proton coincidences from the n + d — > pp reaction. An alternate method used 
elastically scattered pions from 12 C. Pions from this reaction give a strong peak (Fig. [I]) 
with well defined energy for each ring of detectors in the BGO ball. This eliminated the 
strong kinematic energy dependence in the former procedure. Use of the CD 2 as a target 
and a single hit trigger made this possible. The AE and E gains were fitted continuously 
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during data taking and off-line data analysis. This continuous stabilization was important 
because of instability in the BGO ball detectors. The parameters Rij and Z it j also require 
continuous stabilization because both the gain and the decay constant of the BGO scintillator 
is temperature dependent. We used the following stabilization algorithm. For sequential 
subsets of raw events the complete set of parameters was fitted. The weighted average of old 
and new values of the parameters was used for the next subset of events. This procedure was 
completely automated. Fig. || (described below) shows the energy resolution and accuracy 
of this calibration procedure. The width of the two proton total energy peak from the 
d(n + ,p)p reaction, summed over entire ball is a = 3 — 4% MeV. After the system has been 
calibrated one can easily determine the event multiplicity identify detected particles and 
measure their energies. 



IV. RESULTS 
A. Event Selection Criteria 

We selected events using following criteria: 

(1) the multiplicity must be equal 2; 

(2) both particles must be protons; and 

(3) the opening angle between the two protons and their total energy must satisfy the 
kinematics of the reaction n + d — > pp. 

Fig. H] is the energy spectrum for two proton events. The events are summed over all 
kinematically allowed combinations of BGO ball detectors. The narrow peak at the higher 
energy corresponds to the reaction n + d — > pp. The lower energy broad peak with a long tail 
results from pion absorption on carbon. One observes that there is practically no background 
under the n + d — > pp reaction peak. Measurements with no target but with the supporting 
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pipe in its normal position and with a pure carbon target show that the background under 
the peak does not exceed 0.1% so no background subtraction was necessary. 



B. Pion Fraction Determination 

Because of the low beam flux in this experiment it was possible to directly count the 
beam particles. Pulse height and timing information from SI, S2 and DA5 detectors was 
used to determine the pion fraction. First, a cut was placed on the time between SI and S2 
to eliminate accidental coincidences between particles from different beam bursts. Another 
cut was applied to a two dimensional AE — E distribution (pulse height from S2 and DA5) 
to eliminate positrons. The resulting pulse height spectrum from S2 is shown in Fig. |3L The 



left peak corresponds to muons, the right one to pions. Two modified Moyal functions: [11 



F{E) = p^-W-ft)--^-^ (3) 

one for pions and one for muons, were fitted to the energy loss distribution. Here E is the 
energy loss in S2 and Pi — P4 are adjusted parameters. The pion distribution was integrated 
to obtain the pion fraction. At the two lowest energies (21 MeV and 23 MeV) some pions 
loose all their energy in S2 and do not hit DA5 and the corresponding pion peak has a 
different, and more complicated shape. In this case only the muon peak was fitted and the 
pion fraction was determined by subtracting the muon fraction from the integral number 
of beam particles. A Monte Carlo simulation of beam particles passing through the setup 
based on the GEANT code [|12| has been done. The corrections for energy losses, straggling, 
multiple scattering and decay were calculated from this simulation. The resulting correction 
to the total number of beam pions is 7% for the highest beam energy and 44% for the lowest 
one. For beam energies of 21 and 23 MeV this correction accounts for pions, which passed 
through the target but did not hit DA5 and were lost from the beam trigger. This simulation 
also was used for calculation of the pion energy in the center of a target. The uncertainty 
in the number of beam pions, mainly due to fitting errors, is 3 - 5%, with the contribution 
from the statistical error less than 1%. 



C. Efficiency Calculation 



The efficiency was calculated using a Monte Carlo code which incorporates the actual 
geometry of the experimental setup. The efficiency accounts for the effective solid angle 
(geometrical solid angle reduced by the requirement of both protons detection), reaction 
losses, and dead time. The effective solid angles for all kinematically allowed pairs of the 
detectors have been calculated. The missing solid angle due to the nickel cans surrounding 
the crystals is about 6%. Reaction loss correction due to the interaction of protons with 



the BGO were obtained using data from Ref. [|Kj. The typical value of this correction is 
10%. Finally, the data were corrected for dead time. A typical value of the dead time in 
this experiment was 5%. We estimate the total uncertainty in the efficiency of the setup is 
about 5%. 

D. Cross Sections 

The differential cross section has been calculated for the selected events at each specific 
scattering angle. The data for the backward angles were translated to the forward hemi- 
sphere since the angular distribution in the center of mass frame is symmetrical around 90° 
due to indistinguishability of the two protons. The resulting differential cross sections are 
presented in Table |. To obtain the total error, we added the statistical and normalization 
errors in quadrature. An overall systematic uncertainty of about 5% in the efficiency cal- 
culation is not included in the total error. This affects the global normalization of all 90 
differential cross section points and does not change their relative value. 

We fitted the cross sections using a simple parameterization given by the Legendre poly- 
nomial series: 

do 

— = a + a 2 P 2 (cos 6> c . m . ) . (4) 

6^ 'cm. 

At each energy we used only the statistical errors of the cross sections but not the normal- 
ization error because it does not change the relative shape of the angular distribution. 
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The total cross sections were determined by: 

°tot — I "77^ dQ*c.m. (5) 

using the parameterization of Eq. (|]). The integration over 2ti instead of 4n is necessary 
because there are two indistinguishable protons in the final state. The fitted coefficients, 
total cross sections, and y 2 are presented in Table |T|. Errors for ao and 02 are obtained from 
the fit. For the total cross section errors we added the normalization errors in quadrature. 
Overall systematic uncertainty of 5% is not included for the same reason as for differential 
cross section. 

V. DISCUSSION 

The differential cross sections measured in this experiment are plotted in Fig. |j. Solid 
lines show the best fit results with Eq. (|]). The results listed in Table [D] show that only S- 
and P- pionic waves are important in this energy region. The absence of a need for higher 



partial waves agrees with previous measurements |i3| , |T4|| . Fixed angle excitation functions 
are presented in Fig. g. Fig. ^ shows the behavior of the angular distribution as a function of 
energy by plotting the ratio of the coefficients 02 to ao- The shape of the angular distribution 
changes smoothly with pion energy. No structure, within the accuracy of the experiment, 
is observed in the differential cross section as a function of energy. The total cross section 
is presented in Fig. [7| along with the rest of the world's data. The current data follows 
the general trends. The biggest disagreement in the total cross section data, which existed 
around T n = 30 MeV is resolved by the new data. 

We compared our cross sections with the predictions of the recent partial wave analysis 
of Ref. ||15|| . The new VPI solution, SM95, that includes the current data is compared with 



the previous SP94 results which do not include these data. 

Table [TII| and Figs. |5],[7] show the level of agreement between or data and the phase shift 
fits. The normalization factor averaged over all 15 energies are 1.04 and 1.01 for SP94 and 
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SM95 solutions respectively. The values of the normalization factors for both solutions show 
that our estimate of systematic uncertainties is reasonable. Table [IV] presents a comparison 
of our data with other partial wave analysis. The recent VPI results, which include our 
data set, give better agreement than the previous results by the Queen Mary College group 
(9 - 256 MeV) fT§ and by the Hiroshima University group ( 6 - 256 MeV) [|T7J]. The 
good agreement between the predictions of SP94 and recent SM95 VPI solutions with the 
experimental data support the absence of any anomalous behavior of the cross sections in 
this energy region within the accuracy of the present experimental results. 

VI. CONCLUSIONS 

We have presented new data on the total and differential cross section for the reaction 
7T + d — > pp for pion energies from 18 MeV to 44 MeV. In this experiment, 90 experimental 
data points for the differential cross section and 15 points for the total cross section have 
been measured with the same systematic uncertainty. The number of the experimental data 
points in this energy range is almost doubled by these results. The data generally follow 
the trends of previous results and are in a good agreement with the recent SM95 solution 
of the VPI partial wave analysis. The previous disagreement in total cross section data is 
resolved by this experiment. We do not observe any structure either in total cross section 
or in differential cross section. 
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FIGURES 

FIG. 1. Contour plot of the two dimensional distribution of AE (in plastic) vs. E (in BGO) 
obtained from BGO ball phoswich detectors. Regions corresponding to different particle types are 
labeled on the plot. 

FIG. 2. Observed total energy spectrum for two protons produced in pion absorption on a CD2 
target at T w = 43.5 MeV. 

FIG. 3. Pulse height spectrum from the S2 detector for a 30 MeV pion beam. The histogram 
is a measured spectrum, dashed and dotted lines represent the fit results with Eq. (0) for muons 
and pions respectively. Solid line is a sum of these two functions. Positrons have been rejected. 

FIG. 4. Differential cross sections determined in this experiment. Solid lines are the best fits 
with Eq. (g). 

FIG. 5. Energy dependence of the differential cross section at fixed scattering angles. Solid 



lines are the predictions from the recent VPI solution SM95 [15|. 



FIG. 6. Ratio of the coefficients of a Legendre polynomial expansion of the differential cross 
section with Eq. (0j). 

FIG. 7. Total cross section. Black circles are the data from this experiment. Open circles 
represent previous experimental data set for this energy range (taken from SAID |l5| data base. 
Solid line is the VPI SM95 global fit to data in the range - 550 MeV p|. 
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TABLES 

TABLE I. Differential cross sections. Overall systematic uncertainty of 5% is not included in 
the total error. 



T^ ^ era. d(J J d^c.m. 

(MeV) (deg) (mb/sr) 

18.8 36l I 0.891 
41^1 0.793 
58l I 0.669 
66±^ 0.600 
73±ig 0.460 
83±{{ 0.430 

20.9 36l I 0.860 
4ll^ 0.797 
58l I 0.638 
66ljg 0.533 
73t\° 0.437 
83+11 °- 365 

22.9 36l I 0.795 

4ll^ 0.710 

57l I 0.592 

66±$ 0.490 

73t\o 0.456 

84±ii 0.422 

24.0 36l I 0.906 

411^ 0.800 

57l I 0.676 

67±|S 0.621 



statistical normalization total c 

error error (%) error 

0.045 2.9 0.052 

0.035 0.042 

0.038 0.042 

0.033 0.038 

0.030 0.033 

0.028 0.031 

0.043 4.2 0.056 

0.034 0.048 

0.036 0.045 

0.030 0.038 

0.028 0.034 

0.025 0.029 

0.040 4.5 0.054 

0.031 0.045 

0.033 0.043 

0.028 0.036 

0.028 0.035 

0.026 0.032 

0.049 4.3 0.062 

0.038 0.051 

0.040 0.050 

0.036 0.044 
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7311° 


0.458 




84±}i 


0.418 


25.0 


36+1 


0.874 




41+» 


0.808 




58i^ 


0.674 




68±ig 


0.557 




72±ig 


0.443 






0.374 


26.1 




0.984 




411,1 


0.897 




57i^ 


0.728 




67±ig 


0.615 




72±ig 


0.483 




841U 


0.414 


27.1 


361^ 


0.940 




411,1 


0.901 




571^ 


0.757 




6711° 


0.607 




79+10 


0.530 




841H 


0.476 


28.1 


361^ 


0.948 




411,1 


0.931 




571^ 


0.767 




R7+10 

o'-io 


U.DOZ 




7211° 


0.504 




841H 


0.421 


29.2 


361^ 


0.972 



0.032 0.037 

0.029 0.034 

0.045 4.1 0.057 

0.035 0.048 

0.037 0.047 

0.031 0.039 

0.029 0.034 

0.026 0.030 

0.048 3.6 0.060 

0.037 0.049 

0.039 0.047 

0.033 0.040 

0.031 0.035 

0.027 0.031 

0.046 3.7 0.058 

0.036 0.050 

0.039 0.048 

0.032 0.039 

0.031 0.037 

0.029 0.034 

0.049 3.7 0.060 

0.039 0.052 

0.042 0.051 

0.035 0.042 

0.032 0.037 

0.029 0.033 

0.052 3.6 0.062 
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4lti 2 0.965 0.041 0.054 

57l I 0.755 0.043 0.051 

Q7t\o 0.557 0.034 0.039 

72±{° 0.508 0.034 0.038 

Mt\\ 0.512 0.033 0.038 

31.2 36l I 1.032 0.053 3.5 0.064 
4Z±& 0.997 0.042 0.054 
57l I 0.815 0.045 0.053 
67t\° 0.608 0.035 0.041 
72t\o 0.490 0.033 0.037 
85±^ 0.415 0.029 0.033 

33.3 361 I 1.043 0.043 4.2 0.053 
42+^ 0.921 0.032 0.042 
56l I 0.703 0.033 0.039 
67t\o 0.539 0.027 0.031 
72t\o 0.455 0.025 0.029 
85t\l 0.409 0.023 0.026 

35.4 36l I 1.032 0.043 4.0 0.060 
42±& 0.929 0.032 0.049 
56l I 0.768 0.035 0.047 
68t\o 0.611 0.028 0.037 
7lt\o 0.470 0.026 0.032 
85±\\ 0.435 0.024 0.030 

37.4 36l I 1.101 0.068 4.9 0.087 

421^ 0.930 0.048 0.067 

56l I 0.711 0.051 0.062 

68±}2 0.534 0.040 0.048 
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71+^ 0.447 0.038 0.044 

85±{{ 0.408 0.035 0.041 

39.5 36l I 1.205 0.055 3.4 0.068 

42±& 1.171 0.041 0.057 

56l I 0.961 0.045 0.056 

68t\o 0.691 0.035 0.042 

7lt\o 0.567 0.033 0.038 

86±{\ 0.456 0.028 0.032 

43.5 36l I 1.353 0.048 3.8 0.070 

42li| 1.201 0.035 0.057 

561 I 0.946 0.037 0.052 

68±^g 0.732 0.029 0.040 

7lt$ 0.606 0.028 0.036 

86±{{ 0.515 0.025 0.032 

a Energy in the center of target 

b Same value for all six angles at given energy 

c Total error is the statistical and normalization error summed in quadrature. 
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TABLE II. The Legendre polynomial expansion coefficients of the differential cross section as 
fit with Eq. (Q) and the associated total cross sections. Errors for ao and a<2 are obtained from the 
fit. For the total cross section error the normalization errors are added in quadrature. The overall 
systematic uncertainty of 5% is not included. 



T 

- 1 IT 


a Q 


0-2 


X 2 /data 


Otot 


(MeV) 


(mb/sr) 


(mb / sr) 




(mb) 


18.8 


0.663±0.015 


0.458±0.041 


4.7/6 


4.167±0.153 


20.9 


0.632±0.014 


0.510±0.038 


2.4/6 


3.974±0.190 


22.9 


0.597±0.013 


0.372±0.037 


0.7/6 


3.749±0.188 


24.0 


0.669±0.016 


0.478±0.043 


5.5/6 


4.205±0.206 


25.0 


0.649±0.015 


0.517±0.039 


4.1/6 


4.076±0.191 


26.1 


0.717±0.016 


0.583±0.042 


3.8/6 


4.508±0.190 


27.1 


0.730±0.015 


0.503±0.042 


2.0/6 


4.585±0.197 


28.1 


0.732±0.016 


0.576±0.044 


5.7/6 


4.598±0.199 


29.2 


0.743±0.017 


0.557±0.047 


4.1/6 


4.669±0.198 


31.2 


0.763±0.017 


0.683±0.047 


4.9/6 


4.796±0.200 


33.3 


0.714±0.014 


0.652±0.037 


1.8/6 


4.486±0.155 


35.4 


0.740±0.014 


0.624±0.037 


5.9/6 


4.648±0.204 


37.4 


0.724±0.021 


0.690±0.056 


1.5/6 


4.549±0.259 


39.5 


0.888±0.017 


0.834±0.046 


7.3/6 


5.580±0.217 


43.5 


0.939±0.015 


0.857±0.040 


3.2/6 


5.900±0.240 
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TABLE III. The x 2 /data for the present differential cross sections vs pion kinetic energy for 



recent VPI SP94 and SM95 solutions [15|. SP94 does not include the present data while SM95 
includes them. Norm is a common normalization factor determined from the SP94 and SM95 
solutions. 



T n 




SP94 




SM95 


(MeV) 


Norm 


X 2 /data 


Norm 


X 2 /data 


18.8 


1.04 


8/6 


1.03 


7/6 


20.9 


0.98 


3/6 


0.98 


3/6 


22.9 


0.93 


14/6 


0.93 


14/6 


24.0 


1.01 


11/6 


1.01 


11/6 


25.0 


0.98 


8/6 


0.97 


8/6 


26.1 


1.04 


7/6 


1.04 


7/6 


27.1 


1.05 


17/6 


1.05 


16/6 


28.1 


1.04 


12/6 


1.03 


11/6 


29.2 


1.04 


14/6 


1.04 


13/6 


31.2 


1.04 


7/6 


1.04 


7/6 


33.3 


0.96 


4/6 


0.96 


3/6 


35.4 


0.98 


13/6 


0.98 


12/6 


37.4 


0.94 


4/6 


0.94 


4/6 


39.5 


1.08 


13/6 


1.07 


12/6 


43.5 


1.09 


13/6 


1.09 


13/6 


All energies 




147/90 




142/90 
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TABLE IV. Comparison of x 2 /data for the present and World total and differential cross 
sections for the recent VPI SP94 and SM95 solutions [15] and previous solutions BU93 @ and 



HI84 [17]. World data is selected for the energy range 18 - 44 MeV covered by the present 



experiment. 



Solution 


da /dVL c .m. 

Present 


World 


Present 


Otot 


World 


SM95 


142/90 


297/98 


12/15 




38/21 


SP94 


147/90 


292/98 


12/15 




40/21 


BU93 


176/90 


319/98 


52/15 




45/21 


HI84 


332/90 


364/98 


93/15 




73/21 
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